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Abstract 



i-G | It is well known that the matter can significantly alter the picture of 

i ■ neutrino oscillation pj or neutrino decay 0]. Here we show that the pres- 

ence of dense matter induces also the decay of massless majoron, a Gold- 
stone boson associated with the spontaneous lepton number violation, into 
a couple of neutrinos with the same (or in some cases also opposite) he- 
licities. We calculate the rates of such matter induced majoron decays in 
various cases, depending on the neutrino type and the chemical content of 



the medium, and analyse their properties. 
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1. An attractive approach to the neutrino mass problem is based on the con- 
cept of the spontaneous violation of the global B — L symmetry. As a result, 
neutrinos acquire non-zero Majorana masses and the massless Nambu-Goldstone 
appears in the particle spectrum, called the majoron @, [|. It could have inter- 
esting implications for particle physics and cosmology. Since B — Lis free of axial 
anomalies, the majoron is strictly massless and, therefore, is stable. 

As far as B — L is a global symmetry, there can be terms which break it 
explicitly ||. In particular, such terms could appear due to the non-perturbative 
quantum gravity effects, in the form of higher order operators cutoff by the Planck 
scale 0. As a consequence, the majoron x gets non-zero mass, i.e. becomes 
a pseudo-Goldstone boson, and then it could be unstable. Moreover, in order 
to make such a massive majoron cosmologically safe, its decay into neutrinos 
should be necessarily invoked 0. (For the case of the flavon models see also J7|.) 
However, the appearance of the majoron mass is an artefact of some theoretical 
speculations^, and it does not follow from any basic principle. Once again, in the 
framework of the simplest model || the majoron is massless and, therefore, it is 
stable in vacuum. 

In this letter we show that even in the case of strictly massless majoron, its 
decay into a couple of neutrinos (or antineutrinos) x ~ * ^ v ( or X ~ * 2#), becomes 
possible in the presence of enough dense matter background. This phenomenon 
is analogous to the matter induced decay (MID) of neutrino 0. Coherent in- 
teractions with the medium lead to an energy splitting between neutrino and 
antineutrino states, providing available phase space for the emission of the ma- 
joron through the decays v — >z> + x or z/^zz + x, depending on which of v and 
v states becomes "heavier" in the matter background. The neutrino MID and 
its astrophysical implications were subsequently studied in details in a number of 
papers |], [l(J [11], |I2| . However, the fact that the majoron in itself is unstable in 
a medium and can also undergo the MID was missed in these papers.^ 

Here we study in details the dynamics of the majoron matter induced decays 
in various physical situations. The properties of the neutrino field propagating 
in a matter background are given in the appendix. 

2. Let us consider, for the sake of simplicity, the case of one Majorana neutrino 
^> = v + v. The neutrino state is identified with the left-handed spinor v = vl, 
and the antineutrino state v is defined as v = Ur = Cv^. 

The effect of the coherent neutrino scattering off matter constituents can be 
described in terms of the matter potential V. More specifically, we consider a 
neutrino propagating in the rest frame of the unpolarized, constant density matter 
(which is the lab frame in realistic situations). In the framework of electroweak 

^dn the context of certain more complex models, the majoron can acquire non-zero mass 
(without any Planck scale effect) at the quantum level due to axial anomalies. For example, in 
the context of the models B, majoron is simultaneously an axion, in which case its decay into 
neutrinos could be possible. 

2 In fact, such a possibility was first remarked by Rothstein, Babu and Seckel in ref. M. 



standard interactions, the matter potential is given by |I[: 

V = ^pY (1) 

m n 

where Gp is the Fermi constant, m n is the nucleon mass, p is the matter density 
and Y is a number which depends on the neutrino type and the matter chemical 
content (namely, Y = Y e — ^Y n for the u e state and Y = —\Y n for u^ T , where 
Y e , Y n are the electron and neutron number densities per nucleon: Y n = 1 — Y e 
for electrically neutral matter). The potential of the antineutrino state v is —V. 
The neutrino propagation in matter becomes equivalent to the propagation in 
the presence of an external vector field V^ = (V, 0) and can be described by the 
Dirac equation 

{id^ - m + V7V)* = (2) 

Notice, that only the axial neutrino current interacts with the external "matter" 
field Vfj,. Indeed, the the two chiral states in \I/ = ui + vr (which we have inter- 
preted as neutrino and antineutrino) have the opposite " charges" with respect to 
Vfj,. As a result, the energy- momentum dispersion relations are different for the 
positive and negative helicity states (which in the case of the ultra-relativistic 
neutrinos are essentially v and v): 



E + = J(|p|-V0 2 + m 2 (3) 



E- = V(|p| + K) 2 + m 2 (4) 

This modification of the dispersion relations (see e.g. refs. [0]) as compared to 
the case of the propagation in vacuum, allows the 'heavier' of the two helicity 
states to decay into the 'lighter' one with emission of the massless boson.0 The 
natural candidate is the majoron, with the following coupling to neutrinos: 

^x = i 2 X^l 5 ^ = i 2 x("lVr - v R v L ) (5) 

For definiteness, let us consider the case V > (which is e.g. the case of u e 
propagating in normal matter: Y n < 0.66). In this situation the phase space is 
nonzero for the transition v —>■ v + \- The corresponding width is |l~T| , |i~2"| : 



8tt 



W-HH 1 , 



^We < 6 > 



where E is the energy of the initial neutrino. Clearly, when V < 0, as in the 
case of u^ T or in the case of v e in strongly neutronized medium {Y n > 0.66), the 

3 Theeqs. (||), (|]) can be rewritten as E± = |p| 2 +(m^) 2 , where (m^) 2 = m 2 ^2V\i>\+V 2 
are the effective (momentum dependent) 'mass 2 ' terms for the ±1 helicity states. The eventual 
negative rn 2 ,^ (when m 2 < 2V|p|) is not in conflict with the causality principle, since the 
velocity is not proportional to the momentum anymore and it does not exceed the speed of 
light 111. 



v — > z/ + x decay occurs. The neutrino MID exhibites unusual energy dependence 
of neutrino lifetime in lab frame: tmid oc E~ 2 up to E ~ m 2 /V and then, with 
further increasing E, slowly decreases to a constant value T = g 2 V/87r. This is 
just the opposite to the case of usual decay in vacuum where the slow particles 
decay faster: Tvac ^ E. Thus, In the limit of very light neutrino, m 2 <C EV, the 
decay width is determined essentially by the matter potential 0, ||. 

On the other hand, if the matter is sufficiently dense, the majoron decay 
X — > 2z> is also allowed. Indeed, the matrix element is the same as that of the 
neutrino MID, and the phase space considerations are similar. Let us calculate 
now the MID width of the majoron with four-momentum k = (u, k), propagating 
in the medium with constant density. We use the explicit solutions of the Dirac 
equation (^) which are given in Appendix, p = (E', p) and q = (E, q) are the 
four-momenta of the emitted z/'s. We have: 



dT 



22cJ 



\Mf 



(7) 



where the factor 1/2 is due to the identical final states. For the squared matrix 
element, we obtain, neglecting the small terms 0(m 2 /E 2 ) and 0(V/E): 



\M\ 2 = g 2 \C{q,p)\ 2 \a{q) + (3{v)\ 2 - 2 g 2 |p| |q| (1 - cosi 
where 9 is the angle between p and q momenta and 



(8) 



C(q,p) - ^E-\q\ + V)(E'-\p\ + V) 

'(E+\q\-V)(E'+\p\-V) ~ -2V |q| |P| 



(9) 



The differential MID rate is easily obtained by carrying out the trivial integrations 
in the phase space. Due to the four-momentum conservation the r.h.s. of eq. (H) 
becomes 
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so that we have 



dT g 2 V 

dE ~ 8ttu; 
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x 



,uo-E E)\ ' 4Vu 

The physical range of the emitted v energies is 

'l - VI - 4x) < E, E' < - (l + VI -4x 



u> 
~2 



(10) 



(11) 



(12) 



Thus, that the MID x — > 2z> can occur if the majoron energy is above the threshold 
energy Uth = m 2 /V. The total width has the form: 



r 



g 2 V 
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'1- Vl-4x N 
1 + VI - 4x 



(13) 



In the limit i< 1 when the matter potential term dominates over the neutrino 
mass ( i.e. uV 3> m 2 ), the energy distribution of the emitted (anti) neutrinos is 



flat within the energy range ^ < E < oj — ^: 

dT g 2 V 

dE ~ 8nu 

and the total width is independent of the majoron energy: 

t, 9 2 V 

r = - — 

8vr 



(14) 



(15) 



This expression coincides with the one for the neutrino MID v —>■ v + \ in the 
limit of very light neutrinos. For clearness, we remind that the x — > 2z> occurs if 
V > 0. When V < 0, the x ~ y 2^ decay takes place. 

3. Until now we have considered the case of only one Majorana neutrino v 
with matter potential V. Our considerations can be applied to the case of ZKM 
or Dirac neutrinos as well. Let us consider, for example, the case of ZKM neutrino 
^ = u e + v^ resulting from the spontaneous violation of U(l) e (g) 11(1)^ global 
symmetry down to f/(l) e _ At . Here v e = v^ is the electron neutrino (L.H.) state 
with the lepton number L e = 1 and matter potential V e , and v^ = vr = Cz/J 
is a R.H. state of the muon antineutrino with the matter potential — V^. The 
mass term m\I/\I/ = mz/ei/^+h.c. and the neutrino coupling with the majoron 
if X^e^+h.c. conserve the lepton number combination L_ = L e — L^. In the 
case of normal matter we have V e > |V^|. Therefore, the neutrino MID can 
proceed through both the helicity flipping (h.f.) decay u e — ► z> M + \ an d the 
helicity conserving (h.c.) one v e — *■ v^ + \. For the matrix elements of these 
decay modes we have, at the leading order, the following: 

|M| 2 / «2 5 2 |p||q|(l-cos0) (16) 

4 + r7- 2> )(l + cos0) (17) 

q| IpI / 

where p = (E', p) and q = (E, q) are the four-momenta of the initial and final 
fermions respectively and 9 is the angle between p and q. Then the differential 
decay rates are: 

^ + 1 -|) US) 




dT 

dE' 


{v e -> 




dT 



dE' 




(19) 



4 We remind that we have neglected 0(m 2 /E 2 ) terms in \M\ 2 . This is safe since in the 
realistic physical situations, when m ^> V, the final states are always relativistic. If m ~ V, 
one of the final states can be non-relativistic, but the corresponding phase space portion is very 
small (~ V) and it induces only 0(V 2 /E) corrections to the decay width (|i"3|). 



and for the physical range of the final fermion energy we have 

,2 



m 



m 2 + AV±E~ 



-E <E' <E 



(20) 



where V + = \{V e + V^) (for the helicity flipping case) and V_ = \{V e — V^) (for 
the helicity conserving one). The partial widthes are given by: 
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(21) 



(22) 



C J (C + l) 2 . 

Therefore, in the 'small mass' limit (£, (<1) the h.f. decay is dominant, and its 
width depends only on the matter potential: T {y e — > v^x) = 5 ,2 ^ / +/8 7r (cfr- ec L- 



(13)), whereas the h.c. one vanishes as m — > 0. In the opposite regime of 'large 
mass' (£, £ 3> 1) the two widthes have the same shape: 



r(^^) i> e ->!/„*) 2^ 2 e 2 



(23) 



Vf V 3 3 vrm 4 

Both of these are suppressed (by factors £ 3 and ( 3 as compared to the one of eq. 



(13)) and their values are comparable. 

Let us consider the majoron MID channels in the case of ZKM neutrino. 
There is a direct correspondence with the neutrino MID matrix elements due to 
the "crossing" property, li V e + Vn is positive, the majoron MID \ ~ > v e + v„ is 



allowed. The matrix element is still given by the formula 
is given by the formulas similar to (|TT|) and (pT3[): 



and the decay rate 



r(x 



Ve^p) 
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(24) 



Since V^ — V^ > 0, the majoron MID with the opposite helicity final fermions 



is x ~ * v e + Up. Its matrix element is similar to the one of eq. 
differential and total rates are the folowing: 



dT 

dE^ X 



V e Vii) 
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IQ-kuj 2 
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(26) 



Thus, in the limit m — > the h.c. MID rate T(x — > v e v^) is vanishing whereas 
the h.f. one is maximal: T(x — > v e Vy) = g 2 V + /Aii. 

The Dirac case is recovered when the right handed state Vr is sterile i.e. 
Vr = 0. In this case the two majoron MID modes x ~ > v e + Vr and x ~ > v e + vr 
can be directly compared (V+ — V_). 



The straightforward generalization for the case of the majoron matter induced 
decays \ — > V\ + v 2 and \ — > v\ + z/ 2 , where z/ x and z/ 2 are two neutrinos with 
different masses mi,m 2 and different matter potentials Vi, V^, leads to 

r( X -5i^) = ^F(x 1 ,x 2 ), *i=-j^- (27) 

F(xi, a? 2 ) = f{xi, x 2 ) + xxgi(xi,x 2 ) + x 2 g 2 (x 1 , x 2 ) 

n V_ 777, 

r ( X -> ^2) = ^— v^^!, z 2 ) , 2:, = — i- (28) 

47T 4k_w 

~, v /_ mi m 2 \ ,, mi m 2 

G(Zi,Z2) = 2 J{Zl,Z 2 ) gi{Zl,Z2) 92{Zl,Z2) 

V m 2 mi/ m 2 mi 

where 14 = \{V X ± V 2 ) and 

/(a, b) = ^J(l-a-b) 2 -4ab 

/ ,n , fl±(a-b)- f(a,b)\ . . 

4. In order to give some numerical insight on such matter induced decays, 
let us consider the limit of small neutrino mass. Then for the lifetimes of the 
neutrino and majoron MID, we have 

, (l g/cm 3 \ . . 

For example, for the typical supernova densities p ~ 10 14 g/cm , by comparing the 
above equation^ to the diffusive cooling time r ~ Is, we see that these MID effects 
can be relevant down to the values g ~ 10~ 7 , which is well below any laboratory 



and astrophysical limits on the v — \ coupling constants [15], pL6| . However, 
considering the MID effects for the supernova, one has to take into account that 
neutrinos are in equilibrium in the core and comprise a certain amount of the total 
particle number density. Hence they are contributing to the matter potential V in 
eq. (P, as well as providing the Pauli blocking. For example, by taking the case of 
electron neutrino v = u e , we have for Y = Y e — ^Y n + 2Y V = |Yj — ^(1 — Y„), where 
Y[ = Y e + Y v is the density of lepton number trapped in the supernova core. At 
the stage of diffusive cooling, due to the strong neutronization of shocked matter, 
we have typically Y < 0. Therefore the relevant MID's are u e — ► v e + \ and 
X — > 2^ e . However, the formulas @ and (|T3|) cannot be applied directly in this 
case, since one has to take into account the thermal bath of z/ e 's with chemical 
potential \x and temperature T. Therefore, considering either v e — > v e + \ or 
X —> 2i/ e decays, one has to include the factor (1 — f v (E/T)) in the phase space 
per each v e final state, where f v is the Fermi-Dirac distribution function: 

fn\7f) = (E-a\ ; -i ' V = Tf, (31) 

1 exp(^^) + 1 1 



D In fact, in the case of supernova the "small mass" limit £ <C 1 is valid if m < 10 KeV. 



Then for the MID differential rates we have 
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at/ Sttuj 
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(32) 
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(33) 

In both the above equations, the value for the matter potential V refers to the 
'heavier' state in matter, i.e. v e . 

Once again, in the supernova core we can neglect the neutrino mass - m 2 <C 
EV. Then for the MID rates we obtain 



r (z/ e -> u e x) 



8tiE\tJ E Jny T J 



r( x -> 2u e 



g 2 V 



-2 In 

1 



fv^Xe-i + l) 



In 



1 + COsh 7] 

1 + cosh (a) — rj) 



T 




, u = 


UJ 


T 



(34) 



(35) 



8tilu [exp(2?7 — lu) — I] 

In the case of strongly degenerated neutrino thermal background (rj ^> 1), the 
above equations imply for E, Cj < rj 



rfo 



^eX) 



2 V e B - £ - 1 

~8^ W 2 

,2i 



r( x -> 2z/ e 
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whereas for E, uj > r\ we have 



rfe 



^eX) 



r( x -> 2z/ e 
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(36) 
(37) 

(38) 
(39) 
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When i)Cl, and E,u < T, the neutrino and majoron MID are T [y e — > v e \) = 
g 2 V/16ir and T(x — > 2z/ e ) = ^ 2 T^/327r. Hence, the presence of the neutrino thermal 
bath brings just suppression factors 2 and 4 for the neutrino and majoron MIDs 
respectively, as compared to the eq. flip]). When E,u > T, there is no suppression 
and both MID widthes are given by eq. (p~5|). 

5. It is well known that the coherent interaction with matter removes the 
energy degeneracy of different helicity neutrino states. Hence the neutrino, even 
if it is stable in vacuum, can undergo the matter induced decay into the antineu- 
trino and the massless majoron. However, the majoron in itself is unstable in 
matter and decays in a couple of antineutrinos. The matrix element and phase 



space considerations are similar for both the neutrino and majoron MIDs and the 
corresponding widthes are comparable. In this paper the features of the majoron 
MID have been studied for various neutrino types. We have also calculated the 
rates for the neutrino and majoron MIDs in the presence of thermal neutrino 
bath, as it occurs in the supernova core. The impact of the neutrino MID on 



the supernova dynamics and i/-signal was studied in refs. [|, |10|, but the effects 
of the majoron MID were not included. We have seen that these effects can be 
important for v — \ coupling constants g > 10~ 7 -f- 10 -6 . In the 'standard' model 
of the singlet majoron || such a 'large' constant is rather improbable. However, 
there is a variety of singlet majoron models [T7[ which could naturally provide g 



in this range. The effects of the majoron MID for the supernova will be studied 
elsewhere. 

Last but not least, let us remark the following. In fact, the majoron could 
receive considerable mass M, e.g. due to the 'Planck scale' terms explicitly vio- 
lating the lepton number. Then the neutrino MID v — > v + \ could not occur. 
However, in the medium there is still available phase space for the neutrino 3-body 
decay v — > 39 which can be mediated by the exchange of such a heavy 'majoron' 
- or, in other words, due to the lepton number violating 4-fermion interaction 
-ga iyCv) 2 . The properties of such a 3-body MID will be given elsewhere. 



Acknowledgments. 

We thank Georg Raffelt and Dave Seckel for reading the manuscript and useful 
comments. Discussions with Jim Cline, Gianni Fiorentini and Semen Gershtein 
are also gratefully acknowledged. 



Appendix 



We look for the solutions of the Dirac equation (p M 7 M — m + V r 7°7 5 )n(p) = 

in the form of plane waves with negative frequency ^(p) exp(— ipx) and positive 

frequency v (p) exp(+ipx). Let us consider in some detail the former case, where 

we have: 

/ m E + p-a-V \(u\p) \ 

\E-p-a + V m )\ u 2 (p) ) 

Here the Weyl representation for the gamma matrices is adopted; a % are the usual 
2x2 hermitian Pauli matrices. Since the helicity operator ^j is still a constant of 
motion, one can consider the solutions u + (p) and w_(p) with different helicities. 
The presence of the matter potential determines a distinct energy eigenvalue for 
each helicity eigenstate, in contrast with the usual energy degeneracy in vacuum. 



The spinors w±(p) (normalized as u±(p) + u±(p) = 2E±) are the following: 



'E + + |p|-y q(p) \ /-VS_--|p|-V i 9(p) 

m +(p) = h; r~, — 77 , x > u -(p) : 



E + -\p\+Va(p) J \ ^E_ + \p\ + V(3{p) 

where a(p) and /3(p) are the +1 and —1 helicity normalized eigenspinors, respec- 
tively. The matrices a(p)ct(p) + and /9(p)/3(p) + are helicity projection operators: 

a(pWp)+ 4( 1+ lr) ' j9(p)/9(p)+= K 1 ~ipf) 

The normalized positive frequency solutions are obtained in an analogous way: 



, / ^-Ipl + V ^p) \ / _^_ + | p |+y a (p) 

v +(p) = 77; T7 77 „/ s > U ~(P) 



E+ + |p|-^/3( P ) y v-v E --lpi- ya (p 

The energy-momentum dispersion relations for the positive and negative he- 
licity solutions are respectively given by the eqs. (Q) and (f|). The helicity eigen- 
spinors can be chosen in such a way that a(p) and /3(p) are related to each other 
by charge conjugation: 

ia 2 a{pY = (3{p) , ia 2 f3(p)* = -a(p) . 

Summarizing, the (Majorana) solution of eq. @ has the form 

* = *< = ^ ( U+ (p)/ + (p)e-^ +t*_(p)/_(p)«r*-* 
p 

+Mp)/|(p)e +ip+x + ^-(P)/I(p)e +4P - 



where p± = (E±,p) and f± and /|. are the annihilation and creation operators 
(related to each other by charge conjugation) for particle with helicity ±1. There- 
fore, this Majorana field describes the creation and annihilation of two particle 
states only, one characterized by helicity +1 and energy E + (eq. (|3|)), the other 
characterized by helicity —1 and energy E- (eq. (|])). The former is convention- 
ally referred as 'antineutrino' and the latter as 'neutrino'. 
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